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I. INTRODUCTION 



Ferromagnetism and superconductivity are antagonistic to each other states and ferro- 
magnetism, being usually much stronger than superconductivity, destroys the latter. How- 
ever, when dealing with superconductor-ferromagnet {SF) hybrid structures the coexistence 
of these two phenomena is possible. If a normal ferromagnetic metal (F) is in good electric 
contact with a superconductor (S) , the former can acquire some superconducting properties, 
however, with spatially inhomogeneous order parameter [1,2]. The spatial variation of the 
superconducting order parameter in the ferromagnet arises as a response of the Cooper pair 
to the energy difference between the two spin directions in ferromagnet. Due to nonzero 
momentum of the pairs, the usually exponentially decaying Cooper's pair density in the F 
layer has a damped-oscillatory behavior (see, e.g., Ref. [3] for details). As has been theo- 
retically predicted, this causes the formation of the so-called 7r-type superconductivity [4,5], 
an enhancement of dc Josephson current in SFIFS tunnel junctions [6,7] (where / is an 
insulator) and the transition temperature oscillation for SF multilayers [8,9] as a function of 
the exchange field hp or of the thickness dp of a ferromagnet. Experiments that have been 
performed on SFS weak links [10,11] and SIFS tunnel junctions [12] directly prove the 
TT-phase superconductivity. Planar tunneling spectroscopy reveals also a vr-phase shift in the 
order parameter, when superconducting correlations coexist with ferromagnetic order [13]. 
However, experiments on Tddp) SF multilayers behavior are still not so conclusive: the 
nonmonotonic oscillation-like behavior of Tc{dp) was also observed in ,e.g., [14 -16] while 
negative results were reported in, e.g., [17,18]. For interpretation of experimental results, 
along with a mechanism of 7r-type superconductivity and the suppression of Tc due to a 
strong exchange field of a ferromagnet other mechanisms were suggested such as a "magnet- 
ically dead" interfacial layer [16], the effects of a finite interface transparency [18], spin-flip 
scattering [3,19] and other. 

The proximity induced superconductivity of the F metal in SF sandwiches has been 
intensively studied, while not so much attention has been payed to the modification of the 
electron spectrum of the superconductor in a region near the SF interface (an exception is 
a nonequlibrium case, that we do not consider here). In many cases, the Curie tempera- 
ture of the ferromagnet is by order of magnitude larger than the superconducting critical 
temperature Tc , and, in a rough approximation, for SF hybrid structures the "leakage of 
magnetism" into the S layer should be quite strong. Indeed, let us consider the proximity 
effects at the SF bilayer interface. When diffusing into the ferromagnetic layer a supercon- 
ducting electron is subject to an interaction from the local exchange field. E.g., an electron 
with a spin, e.g. "up", has an extra energy ifgxc = A^b^f caused by the intrinsic magnetic 
field hp in ferromagnet {fiB is the Bohr magneton). After tunneling back from the F to S 
layer, one can expect a kind of tunneling interaction that is similar to the superconducting 
proximity effect. Actually, in the superconductor the electron quickly loses its extra energy 
5E ~ ifexc during the time T~h/Hexc being related to a small range Xp~hvp/Hexc is 
Planck's constant, vp is Fermi velocity). Such an equilibrium process leads to a modifica- 
tion of the electron spectrum of the superconductor on nanoscale length Xp . In most cases 
\p~10^^m, so one can assume the importance of such processes for hybrid superconductor- 
ferromagnetic systems with thin S layers. However, in spite of its importance, a theoretical 
analysis of induced magnetic properties of the S metal being in contact with the F one is 
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still lacking. 

The investigation of a " magnetic proximity" effect in SF nanostructures is the purpose of 
our work. On the basis of the microscopic theory for proximity coupled SF bilayer (Sec. II), 
we calculate the spatial variation of the superconducting order parameter into the S layer 
near the interface as a function of the ferromagnetic exchange field and boundary parameters 
(Sec. Ill A). Quasiparticle's densities of states (DOS) variation in the S layer as a function 
of the distance from the SF boundary is also considered on the same basis (Sec. Ill B). It 
is shown that the superconducting pair subjected to the proximity induced exchange field 
breaks, which leads to the formation of quasiparticle subbands in the superconducting energy 
gap. In the absence of spin-orbit (spin-fiip) scattering, the subgap bands accommodate only 
one spin of quasiparticle. Equilibrium local magnetization of the S layer due to the induced 
exchange correlation is also calculated (Sec. Ill C). We end with the Conclusion where we 
emphasize the significance of the leakage of magnetism into superconductor for experiments 
on SF structures. 



Let us consider the proximity effects in the SF bilayer of a massive S {ds ^ ^ Ls) and 
a thin F {Ip dp ^p) layers, for the case of an arbitrary transparency of the SF bound- 
ary, the proximity effect strength and general relation between Curie temperature of the F 



are the thicknesses and Is,f are the electron mean free paths of the S and F layers, re- 
spectively; Hexc is the exchange energy of the ferromagnet. (Henceforth, we have taken the 
system of units with h = ks = 1.) We assume the "dirty" limit for both metals. The su- 
perconducting critical temperature of the F material equals zero. All quantities for S layer 
are assumed to depend only on a single coordinate x normal to the interface surface of the 
materials, while, due to the relation dp <S ^p on thickness, one can neglect the coordinate 
dependence within the F layer. We also expect the F layers magnetization to be aligned 
parallel to the interface, so that they do not create a spontaneous magnetic flux penetrating 
into the S layer. Under these conditions, the only magnetic interaction which can affect 
the superconductor is the short-range exchange interaction between the superconducting 
quasiparticles and magnetic moments into the ferromagnet. 

As is well known, the normal and anomalous Green's functions of the "dirty" metals 
are described by the Usadel equations [20]. We will restrict ourselves to the case when the 
spin-orbit scattering is absent and spin "up" and spin "down" electron subbands do not 
mix with each other (see, e.g., [3]). With the assumptions given above the SF bilayer is 
described by the following equations system for one (e.g." up") spin subband (the domain 
X > is occupied by the S metal and a: < by the F metal) : 



II. PROXIMITY EFFECT MODEL 




$5 = A5 + ,^5' 



(1) 



uGs 



(a;2 + $5$5)i/2' 
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As HT/Tc) + 27rr J][(A5 - ^sGs)/uj] = 0, 

a;>0 



(2) 



= ^T^fCF^p]', Ci.^ . -. (3) 

uGf^ (6^2 + )V2 

Here A5 is the superconducting order parameter of the S material, Cj — uj -\- iHexc, where 
uj = ujn = 7iT{2n + 1), n = ±1, ib2, ±3, ... is Matsubara frequency; the summation over 
frequencies in (2) is cut off by the Debye frequency ljd ; $(a;, i7ea;c) = ^*('^) "-f^exc); the 
prime denotes differentiation with respect to a coordinate x. (Note that the momentum 
renormalization in the ferromagnet is not so important as the frequency renormahzation.) 
We introduce the modified Usadel functions $s = uFs/Gs, $f = loF^^/G^- , where G£,\s 
and Fp^s are Green's functions for the F and S material, respectively, to take into account 
the usual normalized confinement on the Green's functions Gf,s and Ff,s (see Refs. [21,22] 
for details). 

The Eqs. (l)-(3) should be supplemented with the usual boundary conditions in the 
bulk of the S metal: $5(00) = A5(oo) = Ao(T) , where Ao(r) is the BGS value of the 
order parameter, and at the external surface of the F metal $p(— 0?^?) = . The boundary 
conditions at the SF interface are [22]: 



lliGF'<^'F 
UJ 



x=0 



LO 



(4) 



{7bfGf^>fL=o = ^Gs{^s/(^ - ^F/i^)\^=Q , (5) 

where 7 = Ps^s/Pr^ is a measure of the proximity efi^ect strength, '-/bf = Rb/Pf^ describes 
the effect of the boundary transparency. Here ps,F are the normal-state resistivities of the 
S and F metals, Rb is the product of the SF boundary resistance and its area. Here 

and below we write our formulas for the F metal using the effective coherence length of 

/ \ 1/2 

normal nonmagnetic (N) metal C = [ 2:^^ ) with the diffusion coefficient Df instead of 

/ \ 1/2 

~ ( ) ' have a possibility to analyze both limits Hexc {SN bilayer with 
(In <^ ^) and H^xc >> ■ The relations (l)-(5) generalize the proximity effect problem 
with an arbitrary interface transparency for the case of normal metal with ferromagnetic 
order. 



III. RESULTS 



A. The spatial variation of the order parameter 

Due to a mesoscopic thickness of the F metal the proximity effect problem can be reduced 
to the boundary problem for the S layer and a relation for determining [22,23]. There 
are three parameters which enter the model: 'Jm — idp/^ is the measure of the strength 
of proximity effect between the 5* and F metals, 7^ = iBFdF/i describes the electrical 
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quality of the S/F boundary, and H^xc is the energy of the exchange correlation in the F 
layer. In the general case, the problem needs self-consistent numerical solution of the Usadel 
equations (such mode consideration see, e.g., in Ref. [24]). Here, however, we will not discuss 
the quantitative calculations, but present qualitatively correct ones to consider new physics 
we are interested in. In this approximation, for bilayer with weak proximity effect 7m <S 1 
and arbitrary value of the boundary resistance 7^ (quite a realistic experimental case) for 
the function $5 we have the analytical solution (see Ref. [22,23]): 

M^..) = A„{l-,«;3^ ^^<-f/fe\ } (6) 

^Mpoj + ujA[uj) 

where = {u^ + l^lf'^ /ttTc, and A(cu) = [1 + ^bCj {ibCj + 2uj / p^) /{TiTcff'\ Using 
$5(0;, a;) (6) and the self-consistency condition (2) one can obtain the spatial variation of 
the order parameter in the S layer /^s{x) for different values of 75 and 7m <^ 1. As follows 
from the analysis, the BCS value of the order parameter is reached at the distance of several 

into the S layer. At fixed value of the exchange interaction Hexc , the order parameter 
at SF boundary As{0, Hf^xc) decreases with increasing boundary transparency, which leads 
to a decrease in the jump of the amplitude of the Cooper pairs in going from the S to the 
F layer. A decrease in the boundary transparency leads to a sharp increase of the order 
parameter at SF boundary. The behavior on 7m and 75 parameters here is analogous to 
the situation in the SN system (see, e.g., Ref. [21]). 

Another important feature of Eq.(6) is its dependence upon Hexc ■ The exchange inter- 
action influence on the spatial variation of the order parameter in the S layer is shown on 
Fig. 1. Namely, the curves show the dependence of difference of the order parameters in the 
S layer for the case, when magnetic interaction is turned off {SN bilayer) and the case with 
ferromagnetic correlation {SF bilayer) on distance from interface for fixed boundary parame- 
ters and different strength of the F layer exchange field H^xc '■ [^s{x, H^xc = 0)— A5(a;, Hexc)] 
. It is seen that magnetic influence decreases with the increasing of the distance from SF 
boundary and at x > 2^5 the order parameter becomes equal to those for the case ol SN 
bilayer. This is not surprising, since one would expect that the induced exchange field 
should suppress the superconducting order parameter at some distance into the S layer in 
excess of the layer suppression. Suppression increases with the increase of the exchange 
interaction parameter ifea;c and of the proximity effect. The curves on Fig. 1 illustrate the 
spatial dependencies of the induced exchange correlation in the superconductor for the case 
of vanishing interface resistance 75 = . With an increase of the SF boundary resistance 
the electrical coupling of the S and F metals decreases and in the limit 75 — 00 the 
metals become decoupled. As follows from (6), far {x ^ ) from the interface the bulk 
superconductivity is re-established. 



B. Spin-splitting of DOS and local exchange-interaction induced states 

Spin sphtting of DOS in the S layer is the second demonstration of a magnetism leakage 

into the superconductor. The Green's functions for the S layer Gs]]{i^^x) and Gsii{i^^x) 
for both spin subbands one can obtain using the solution (6) for the functions ^^(cj, x) with 
Hi = ui + iHexc and lj = lu — iH^xc , respectively. Performing the analytical continuation 
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to the complex plane by the substitution uj —> —is wc calculate the spatial dependence 
of quasiparticle DOS for each spin subband Ns]{s,x) = ReGs]]{i^.x) and Nsi{£^x) = 
ReGsii{u!,x). On Fig. 2 the Ns^{e,x) dependence at different distances from the SF 
interface is presented for T << T^, i^eic = SttT^ and jm = 0.1, and vanishing boundary 
resistance (7^ = 0) . We find that for the case of Hexc 7^ 0, and 7m 7^ the quasiparticle 
DOS is spin-splitted, that means Ns^{s,x) 7^ Nsi{e,x) and Ns^{e,x) ^ Nsi{—s,x), in the 
S layer at the distance of a few C,s from the SF boundary. As one can expect from the 
fermionic symmetry, the spin-up particles and spin-down holes have the same DOS, and 
likewise for spin-down particles and spin-up holes. The spin-splitting decreases with an 
increase of the distance from the boundary and vanishes in the bulk of the ^S" layer (see curve 
4 on Fig. 2). The DOS of Cooper's pairs = ReFs^ 1(00 , x) , Kj^(e,a;) = ReFsi^{uj,x) is 

spin-splitted, too. These features are due to the initial spin-splitting electrons on the Fermi 
surface in the F metal and characterize the SF bilayer as a united system. 

On Fig. 3 the Ns'^ie^x) dependence is presented for = 1 and different values of 
the exchange energy for bilayer with finite boundary transparency 7^ = 0.1 and proximity 
effect strength jm = 0.1. The behavior on 7^ parameter is due to two mechanisms: with 
decreasing boundary transparency spin-splitting increases if 75 ^ 0.2, however, for large 
enough ■jb > 0.2 the decrease of the interface electrical quality preponderates and DOS 
approaches the BCS value. 

Other important features, shown on Figs. 2 and 3, are the local states that exist inside 
the energy gap at the distances of a few from the SF boundary. These subgap states are 
absent at the SF boundary {x = 0), far from the SF interface, and also absent if Hexc = 0. 
For small values of jm and 75 = , as follows from Eq.(2) and solution (6), Ns-i{s,x) has 
singularities for 

e = ±Ao{l - l^'^'l . eM-Psx/^s)} (7) 

where (3^ — (Aq — e^Y/'^ /t^Tc and e = e — Hexc- We found the singularities inside the 
superconducting gap, — Aq < e < Aq , by numerical calculations. The solutions of Eq.7 
for x/^s = 0, 1,5 and for 7^ = 0.1, Hexc = ^t^Tq are presented on Fig. 4 in the diagram 
form. The local states are thus definitely not due to the spatial variation of the pair po- 
tential, but due to Cooper pairs breaking in the superconductor by the exchange-induced 
magnetic correlation. The region of their existence increases with the increasing of Hexci or 
increasing pair breaking effects. In the absence of spin-ffip (e.g., spin-orbit) scattering, the 
subgap bands accommodate only one spin of quasiparticles. These bands bear superficial 
resemblance to both the bands observed at interface of superconductor and perfectly insulat- 
ing ferromagnet [25] and bulk superconductor containing finite concentrations of magnetic 
impurities [26,27]. 

C. Local magnetization 

As was mentioned above, the infiuence of the ferromagnet on the superconductor is re- 
fiected in a nonzero value of the difference in the densities of states for spin-up and spin-down 
electrons, Ns^ {e, x) and Nsi {e, x) . This DOS splitting causes inhomogeneous quasiparticle 
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spin density in the S metal, i.e., an effective magnetization Ms{x) of the S layer, that can 
be found using the relation (see, e.g., [28]) : 



where Mq = gSeiisi^ IJ^b) , -S'e = 1/2,6/ = 2 and f{e) = l/[exp(£/r) + 1] is Fermi 
distribution function. This is confirmed by numerical calculations of Ms{x) Eq.(8) shown in 
Figs. 5-7. Fig. 5 shows the leakage of magnetization into the superconductor versus distance 
from the SF interface for fixed boundary parameters and different exchange field energy 
values. Fig. 6 shows the same magnetic characteristics but for SF sandwich with fixed 
exchange energy and boundary transparency, and different proximity effect strength. The 
curves on Fig. 7 illustrate inhomogeneous magnetism induced into the superconductor by 
proximity effect, Ms{x) , for the SF bilayers with different interface transparency and fixed 
i^eic and 7m • Following DOS dependance, the magnetization behavior on 7^ parameter 
is due to two mechanisms: firstly, 7^ ^ 0.2 , with decreasing boundary transparency spin- 
splitting increases and Ms{x) increases too, secondly however, for large enough 7^ > 0.2 
the decreasing of the interface electrical quality preponderates and induced magnetization 
decreases. 

All the curves of Ms[x) show essentially non-linear behavior versus distance from the SF 
interface, and have important common feature - a damped oscillatory behavior on a scale of 
superconducting coherence length ^5 . As is evident from the plots on Figs. 5-7 , induced by 
proximity effect magnetic correlations in the S layer are inhomogeneous in character with 
spatially dependent magnitude and direction of its molecular field. 



We have calculated the magnetic proximity effect induced in the superconductor at SF 
interface. We find that an equilibrium exchange of electrons between the F and S metals 
results not only in proximity induced inhomogeneous superconductivity of the F metal, as 
was found earlier, but in proximity induced inhomogeneous magnetism of the S metal, too. 
Our results show that the appearing damping magnetic correlations in the superconductor 
are quite robust and, to our opinion, are quite possible for experimental measurements. 

Experimental studies of magnetic proximity effect in SF structures are nowadays unclear 
and controversial. Tedrow et. al. [29] were the first who carried out high-field tunneling 
measurements of the quasiparticle DOS in superconducting films of Al backed by EuO, 
which in bulk form is a ferromagnetic insulator with a Curie temperature Tcr = 70K. The 
tunnel conductances of these junctions show a spin splitting of the Al DOS in excess of 
the splitting observed for Al films in contacts with nonmagnetic substrates in equivalent 
magnetic field. The authors interpreted the excess spin splitting in terms of a coupling 
between the Al conduction electrons and the EuO magnetization. However, Tokuyasu et.al. 
[30] argued that the observed excess Zeeman splitting in EuO /Al proximity contacts cannot 
be explained in such a mode, and considered that excess splitting is the reflection of the 
exchange coupling of the tunnehng electrons with Eu^'^ spins at the interface. Recent 
experiments [31,15] revealed that the effective thickness of the magnetic layer in SF hybrid 
structures usually is much larger than its physical thickness. The increase of the thicknesses 
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IV. CONCLUSION 
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is so great that in all samples, except for those with extremely thin magnetic layers, the 
crossover to a 3D state superconductivity is never in fact observed experimentally. This is 
to be contrasted with the case of nonmagnetic spacer layers, where these two length scales 
are comparable. Taking into account our results, we explain that the rise of the effective 
magnetic layer thickness in the SF sandwiches is due to magnetic proximity effect. Namely, 
the induced magnetic correlation into the S layer depletes Cooper's pair density at the SF 
boundary that results in an excess thickness of the magnetic layer. 

In conclusion, we study the magnetic correlations into a superconductor that have been 
induced at SF interface due to proximity effect. Equilibrium leakage of magnetism into the 
S metal results in (i) additional spatial suppression of the order parameter in the S layer, 
(ii) spatial spin splitting of the quasipaticles DOS, (iii) local bands that appear inside the 
energy gap in the S layer, and (iv) spatially dependent magnitude and direction induced 
magnetization of the S layer. The existence of these magnetic properties of the S metal is 
quite important for SF nanostructures and should be taken into account while comparing 
theoretical results with experimental data. 
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Figure Captures 

Fig. 1. The difference of the superconducting order parameter in the S layer versus 
distance from the interface for SN and SF structures with the same boundary parameters 
(7m = 0.1 , 7b = ), and different ferromagnetic field energy H^xc/t^Tc =8, 9, 10, 12 and 
15. 

Fig. 2. Normalized density of state for spin "up" quasiparticles in the S layer of the 
SF sandwich for 7jv/ = 0.1 , 7^ = and iJe^rc = SttT^ , and various distances from the SF 
interface: x/^s = Oi li 5, and 30 (curves 1, 2, 3, and 4, respectively). 

Fig. 3. Same as on Fig. 2 for 7^/ = 0.1 , 7b = 0.1 and x = ^5 , and various ferromagnetic 
field energies: Hexc/'^Tc =1, 2, and 5 (curves 1, 2, and 3, respectively). 

Fig. 4. Diagram solutions of the Eq.7 for 7m = 0.1, H^xc — ^t^Tc and x/^s — 0, 1, 5, 
(curves 1, 2, and 3, respectively). The regions where the local states exist are shaded; see 
also Fig. 2 for comparison. 

Fig. 5. Leakage of magnetization into the S material versus distance from the interface 
for SF sandwich for 7m = 0.1 , 75 = , and different exchange energies H^xc/t^Tc =7, 5, 
and 3 (curves 1, 2, and 3, respectively). 

Fig. 6. Same as on Fig. 4 for 75 = , Hexc — 3.57rTc and different proximity effect 
strength 7m = 0.1, 0.15, 0.2 (curves 1, 2, and 3, respectively). 

Fig. 7. Same as on Fig. 4 for 7m = 0.1 , i/exc = SttTc and different SF boundary 
transparency 7^ = 0.0, 0.05, and 0.1 (curves 1, 2, and 3, respectively). 
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